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Abstract 

In the general case of the so(n) free rigid body we will give a list of integrals of motion, which 
generate the set of Mishchenko's integrals. In the case of so(5) we prove that there are fifteen 
(-H ' coordinate type Cartan subalgebras which on a regular adjoint orbit give fifteen Weyl group orbits 

^ (— I of equilibria. These coordinate type Cartan subalgebras are the analogues of the three axes of 

equilibria for the classical rigid body on so (3). The nonlinear stability and instability of these 
equilibria is analyzed. In addition to these equilibria there are ten other continuous families of 
equilibria. 
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1 Introduction 

qq The purpose of this paper is to study the geometry and the dynamics of the free rigid body on the Lie 

CO ■ algebra so (5). The general free rigid body on so(n) has been studied as a completely integrable system 

in the classical works [T], [2J, [3]. The general case and more extensively the even case n — 2m have 
Q\ ■ been analyzed by Feher and Marshall [3] , Spiegler [5] , where a method for determining a certain class of 

equilibria is given. They also study the stability of these equilibria using the Energy-Casimir method. 
As mentioned in [5], the odd case n — 2m + 1 is significantly different and only some indications on 
stability results have been given. 

fn [6J it has been proved that this class of equilibria comes as the intersection of coordinate type 
Cartan subalgebras with regular adjoint orbits. Consequently these equilibria are unions of Weyl orbits 
and they correspond to long-short axis type of equilibria known from the dynamics on so (3). 

The stability of these equilibria, for the case of so (4), has been studied in [5J, using Williamson normal 
form [TJ. For the case of so (4) there have been discovered two new continuous families of equilibria on 
every regular adjoint orbit. Also, in <6 it has been shown that these are nonlinearly stable as a family, 
that is, if a solution of the so(4)-free rigid body equation starts near an equilibrium on such a curve, at 
any later time it will stay close to this curve but in the direction of the curve itself it may drift. 

First of all, in this paper we will give a list of integrals of motion for the general case of the so(n) free 
rigid body, which proves to be a set of generating functions for the Mishchenko's quadratic integrals of 
motion and which will play a crucial role in the approach of the stability problem using energy methods. 

For the case of so (5) rigid body we found fifteen families of coordinate type Cartan subalgebras, 
which intersected with a regular adjoint orbit give 120=15x8 equilibria corresponding to long-short 
axis type of equilibria (8 being the cardinal of the Weyl group of so(5)). We also found ten continuous 
families of equilibria on every regular adjoint orbit. As in the case of so (4), these two types of equilibria 
give all the equilibria of so (5) free rigid body. 

The stability problem is studied using the linearization method and Arnold's method, which is 
equivalent to Energy-Casimir method [S]. In order to do this, we explicitly compute the integrals 
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of motion using Mishchcnko's and Manakov's methods and we point out the need for the additional 
constants of motion introduced in Section 2 of this paper. Extended indications on how Arnold's method 
applies for each studied equilibria are given. 



2 A new family of integrals of motion for the free rigid body 
on so(n) 

The equations of the rigid body on so(n) are given by 

M = [M, Q] , (2.1) 

where ft £ so(n), M — VLJ + Jil G so(n) with J = diag(A,), a real constant diagonal matrix satisfying 
Aj + Xj > 0, for all i, j = l,...,n, i ^ j (see, for example, [3]). Note that M = [rriij] and 57 = [wy] 
determine each other if and only if Aj + Aj > since rriij = (A, + Xj)u>ij which physically means that 
the rigid body is not concentrated on a lower dimensional subspace of R" . 

It is well known and easy to verify that equations (|2.ip are Hamiltonian relative to the minus Lie- 
Poisson bracket 

{F, G}(M) := - Trace(M [VF(M), VG(M)]), (2.2) 
and the Hamiltonian function 

H(M) := --Trace(MQ). (2.3) 
Here F,G,H S C°°(so(n)) and the gradient is taken relative to the Ad-invariant inner product 

(A, Y) := —- Trace(Ar), X, Y e So(n) (2.4) 

which identifies (so(n))* with so(n). This means that F = {F,H} for all F G C°°(so(ri)), where 
{•,•} is given by (|2.2I) and 7? by (I2.3[) . if and only if (|2.ip holds. Note that the linear isomorphism 
A G so(n) h-> XJ+JX G so(n) is self-adjoint relative to the inner product (|2.4p and thus VH(M) = O. 

In what follows, we assume that are all distinct. As stated in [T], the system (I2.ip admits a 
sequence m r (r = 1, 2, 3, . . . ) of integrals of motion, all of them depending quadratically on the angular 
momentum M. It is easy to compute an explicit form for these integrals [9] as follows: 

m r (M) = J2 ^Arn^ k . 

i<k 

Next we will find an additional family of n integrals of motion for the system (|2.1I) . which have simple 
and elegant expressions and which prove to generate Mishchenko's integrals of motion. The complete 
result is contained in the following 

Theorem 2.1. The functions 

1 

Fi{M) = J2 A2 _ A2 m 2 ik , i = l,2,...,n 

fe=i i k 

are integrals of motion for the system (12. ip ; moreover, the following identities hold for each r: 



i r (M)=J2>H F i(M). 
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Proof. It is easy to see that 



•(ao = E 



i,fc=l 
i<k 



E 



mj 



Af-Ar 



'ife 



E 

i,fc = l 
n 

E 

i,k = l 

k<i 



A 2 — A 2 



Af-Ar 



n _^ 71 1 

fc=l 4 K fc=l ' K 



\k>i 



k<i 



n n \ 

i—1 k=i * ^ 

k^i 

n 

= E A ^w- 



Now, writing the above identities for r = 1, 2, . . . , n we obtain a Vandermonde-type linear system, which 
solved leads to the conclusion that the functions Fi are linear combinations of Mishchenko's integrals 
of motion m r (r = 1, 2, . . . , n) and thus are integrals of motion for our system (12.11) . ■ 



3 The free rigid body on the Lie algebra so (5) 

We begin by making some considerations on so(5), which is the Lie algebra of the subgroup SO(5) 
{A e g[(5,R) | A t A = 7 5 ,det(A) = 1} of the special linear Lie group SL(5,R). 
We choose as basis of so (5) the matrices 
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and hence we represent so (5) as 
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Relative to this basis, the Lie algebra structure of so (5) is given by the following table 
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In the basis {Ei, . . . , Eiq}, the matrix of the Lie 
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structure (|2.2[) is 
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Since rankso(5) = 2, there are two functionally independent Casimir functions which are given respec- 
tively by 

Ci(M) := -\ Trace(M 2 ) = \ ( £ x\ + £ + £ 

\i=l j=l i=l / 

and 

C 2 (M) := iTrace(M 4 ) = 
8 

= 5 [(^2 + x\ + y\ + z 2 ) 2 + (x 2 + x 2 + y 2 + z 2 ) 2 + (x 2 + x 2 + y\ + z 2 ) 2 + 

+(2/i + 2/2 + 2/3 + zf) 2 + (^i 2 + zl + 2| + zf ) 2 + 2(yizi + 27 2 z 2 + y 3 z 3 ) 2 + 
-\-2{x 2 y 3 - 2:32/2 - 2az 4 ) 2 + 2(x 3 j/i - xiy 3 - z 2 z 4 ) 2 + 2(xiy 2 - x 2 yi - z 3 z 4 ) 2 + 
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+2(x 3 z 2 - x 2 z 3 - yiz 4 ) 2 + 2(x x z 3 - x z z x - y 2 z 4 ) 2 + 2{x 2 z 1 - x x z 2 - y 3 z 4 ) 2 + 
+2(x 1 x 2 - yiy 2 - z x z 2 ) 2 + 2(x 1 x 3 - y x y 3 - z 1 z 3 ) 2 + 2(x 2 x 3 - y 2 y 3 - z 2 z 3 ) 2 ] . 
Thus the generic adjoint orbits are the level sets 

Orb ClC2 (M) = (Ci x C 2 )-\c 1 ,c 2 ), (ci,c 2 ) G K 2 . 

In all that follows we will denote by Orb Ci;C2 the regular adjoint orbit Orb ClC2 , where c\ > 0,c 2 > 
and 2c 2 > c\> c 2 . 

Using the Lie bracket table in the chosen basis given above, it is immediately seen that the coordinate 
type Cartan subalgebras of so(5) are tk, 1 < k < 15, where 



ti := span(E 3 ,E 6 ) 



<6 



t 2 := span(E 6 ,E$) = < 



t 3 := span(£; 6 ,£;7) 
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t 4 := span(E 5 ,^ 7 ) = < 



<6 



t s := span(£i,£V) = < 



M. 
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ie := sp&n(E 2 ,E 5 ) = < 



t 7 := span(£; 5 ,£;9) = < 



<6 



t 8 := span(^i,Si ) = < 
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t g := span(£J 1 , E 4 ) = < 



ML : 



tio := spa,n(E 2 ,E 8 ) = < 



tn := span(S 4 ,£; 8 ) = < 



ti2 := spa,n(E 3 ,E w ) = < 



M. 
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ti 3 := span^,^) = < 
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ti4 := span(i?4, _Eg) = < 
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The intersection of a regular adjoint orbit and a coordinate type Cartan subalgebra has eight elements 
which represents a Weyl group orbit. Specifically, we have the following result: 



Theorem 3.1. The following equalities hold: 

ik n Orb Ci;C2 = {Ml b , M\_ h , Mt a , M\_ a , M\ h , M^ bl 
for all k, 1 < k < 15, where 



,M\ a } 



a = \/ci + ^ 2c 2 - c\ 
b=\jd- \ 2c 2 - c\. 



(3.3) 
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Proof. Let Etjfl Orb Ci;C2 . Then € h, c x = C^M^) = 1 (a 2 + (3 2 ) and c 2 = C 2 (M^) = 

j (a 4 + (3 4 )- We obtain the system 

/ a 2 + /3 2 = 2 Cl 
\ a 4 + (3 4 = 4c 2 

which leads immediately to the result. For ik n Orb Ci;C2 , 2 < k < 15, we proceed in a similar manner. ■ 



4 Equilibria for the so(5)-rigid body 

We will work from now on with a generic so(5)-rigid body, that is, Xi + Xj > for i =/= j and all Aj are 
distinct. The relationship between ft = 6 so(5) and the matrix M 6 so(5) in the representation 
(|3.ip is hence given by 



(A 3 + A 2 )w32 = xr, (Ai + A 3 )o;i3 = a; 2 ; (A 2 + Ai)w 2 i = ^35 

(Ai + = yr, (A 2 + A 4 )w 24 = ?/2 ; (A 3 + A 4 )w 34 = y 3 ; 

(A x + A 5 )u; 15 = zx; (A 2 + A 5 )w 25 = z 2 ; (A 3 + A 5 )w 35 = z 3 ; (A 4 + A 5 )w 45 = z 4 . 
and thus the equations of motion (|2.1I) are equivalent to the system 



xi = (A 2 - A 3 ) 

±2 = (A 3 - Ai) 
X3 = (Ai - A 2 ) 
2/i = (Ai - A 4 ) 
2/2 = (A 2 - A 4 ) 
2/3 = (A 3 - A 4 ) 
zi = (Ai - A 5 ) 
z 2 = (A 2 - A 5 ) 
Z3 = (A3 - A 5 ) 
i 4 = (A 4 - A 5 ) 



112113 



(A 2 +A 4 )(A 3 + A 4 ) (A 1 +A 2 )(A 1 + A 3 ) (A 2 + A 5 )(A3 + A 5 ) 



(Ai+A 4 )(A 3 +A 4 ) (Ai+A 2 )(A 2 + A 3 ) 1 (Ai+A 5 )(A 3 +A 6 ) 



V1V2 



+ 



(Ai+A 4 )(A 2 +A 4 ) (A 2 +A 3 )(A! + A 3 ) 1 (Ai+A 6 )(A 2 +A 5 ) 



X2V3 



(A 1 +A 3 )(A 3 +A 4 ) (Ai+A 2 )(A 2 +A 4 ) (A!+A 5 )(A 4 +A 5 ) 



(Ai+A 2 )(A!+A 4 ) (A 2 +A 3 )(A 3 +A 4 ) (A 2 +A 5 )(A 4 +A 5 ) 



■ run 
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(A 2 +A 3 )(A 2 +A 4 ) (Ai+A 3 )(Ai+A 4 ) (A 3 +A 5 )(A 4 +A 5 ) 
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(A1 + A3KA3+A5) (Ai+A 2 )(A 2 +A 5 ) ^ (A!+A 4 )(A 4 +A 5 ) 



(A 1 + A 2 )(A 1 +A 5 ) (A 2 +A 3 )(A 3 +A 5 ) ^ (A 2 +A 4 )(A 4 +A 5 ) 



K3Z4 



(A 2 + A 3 )(A 2 +A 5 ) (Ai+A 3 )(Ai+A 5 ) (A 3 +A 4 )(A 4 +A 5 ) 



!l2~2 



(A 1 +A 4 )(A 1 +A 5 ) (A 2 + A 4 )(A 2 + A 5 ) (A 3 +A 4 )(A 3 + A 5 ) 

The Hamiltonian (|2.3p has in this case the expression 
1 



H{M) = --Trace(Mft) 
1/1 



2 V A 2 + A 3 1 Ai + A 3 2 Ai + A 2 3 Ai + A 4 W1 A 2 + A 4 



-vi + 



A 3 + A 



■2/3" 
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1 



Ai + A 5 A 2 + A 5 A 3 + A 5 A 4 + A 



(4.1) 
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The Hamiltonian nature of system (|4.ip can be checked in this case directly, writing 

(±1, x 2 ,± 3 , yi,y2, i/3,h, z 2 ,z 3 , i 4 ) T = r_(Vi?) T , 

where the Poisson structure T_ is given by (|3.2[) . 



/ 15 \ / 10 \ 

Theorem 4.1. If £ denotes the set of the equilibrium points of (|4.1I) . i/ien £ = I [J 1 U I (J S; ), 

where Si, 1 < Z < 10, are £/ie i/iree dimensional vector subspaces given by 

Si, 2 := span R { (j^I^i ± x£% E *) , [ T ^E 2 ± x^^s) , ( XITaI^ ± aTTaI^) } ; 

s 3 ,4 := span K { (j^Ei ± i^^io) , ( a^Fa^s ± a^h^) , ( jjfjj E e T a^aI^s) } ; 

s 5 , 6 := span R { (x^£ 2 ± j^; E io) , ( aTTa^g ± j^ e t) , ( j^M E 4 T aTTaI^s) } i 

s 7 ,8 := span R { ( j^^i ± a^FaT^) > ( a^Fa^ 2 ± x^xl E s) , ( aITaI^s ± aITa;^) } i 

s 9 ,io := span R { (j^t^s ± aTTa^io) > (a^^ 4 ± aTTaI^s) - ( aTTaT^s T a^FXI^) } . 
Proof. The proof follows after a long, but straightforward computation. ■ 

5 Constants of motion and nonlinear stability 

In this section we study the nonlinear stability of the equilibrium states £ f~l Orb Ci;C2 for the dynamics 
(14. ip on a generic adjoint orbit. 

Using Mishchenko's method [T], [3J, [5], we obtain the following additional independent constants of 
the motion for Eqs. (|4.1[) : 



1 3 
Jfi(M) = — Traced J p MJ 3 - p n 



4 

p=0 

- i[(Al + A§)a^ + (A? + A 2 ):r| + (A? + \%)x\ + (A 2 + + {X% + X\) v \ + (A3 + A 2 ) % 2 
+ (A? + A 2 )z 2 + (A 2 + xl)zl + (A§ + A 2 )z 2 + (A 2 + A 2 )zf] 

and respectively 

1 5 

# 2 (M) = — Trace ^ J p MJ 5 ~ p n 

p=0 

= \m + 4 + 44)4 + (4 + 4 + 44)4 + (4 + A 4 + A 2 A 2 )x 2 + 

+ (A 4 + A* + A 2 A 2 )?/ 2 + (A 4 + A 4 + X 2 2 X 2 4 )y 2 2 + (A 4 + A 4 + A 2 A 2 ) % 2 + 

+ (A 4 + A 4 + A 2 A 2 )z 2 + (A 4 + A 4 + A 2 A 2 )z 2 + (A3 + A 4 + A 2 A 2 )z 2 + (A 4 + A 4 + A 2 A 2 )zf], 

which are, more precisely, the Mishchenko's integrals of order 4 and respectively 6. 
Thus Mishchenko's method provides two constants of motion (the above constants K\ and K 2 ), which 
add to the Hamiltonian H. As mentioned in [3J, the number of independent constants of motion 
generated by this method is half of the dimension of the adjoint orbit only for so (4). It follows that 
using Mishchenko's method we cannot obtain other independent constants of motion for the rigid body 
on so (5) except the already mentioned ones. 



In what follows, we will analyze the results of Manakov's method for finding constants of motion ([2], 
[3], [5]), which are to be found as coefficients of the powers of 7 in the expansion of 

— TracefM + 7 J 2 Y 
2r 

for r = 2,3,4,5. 

For r — 2 we obtain the first Casimir, C\. For r = 3we obtain the constant of motion K\. For r = 4we 
obtain the constant of motion K 2 and the second Casimir, C%. Finally, for r — 5 we obtain a constant 
of motion which can also be obtained by Mishchenko's method for order 8, but is dependent of the 
previously obtained ones, and the new functionally independent constant of motion 

where 

T\ = {x 3 z 2 - x 2 z 3 - yiz 4 ) 2 + (x 2 y 3 - x 3 y 2 - z x z 4 ) 2 + (x x x 3 - y x y 3 - z 1 z 3 ) 2 + 

1 / \2 1 / 2 1 2 1 2 1 2\2 

+ (xix 2 - yiy 2 - Z1Z2) + (x 2 + x 3 + y 1 + z x ) ; 

T 2 = (xiz 3 - x 3 zi - y 2 z 4 ) 2 + (x 3 yi - x\y 3 - z 2 z 4 ) 2 + (x 2 x 3 - y 2 y 3 - z 2 z 3 ) 2 + 
+ (xix 2 - y x y 2 - z 1 z 2 ) 2 + [x\ + x 2 3 + y 2 + z 2 ) 2 ; 

T 3 = {x 2 zi - x x z 2 - y 3 z 4 ) 2 + {xiy 2 - x 2 y x - z 3 z 4 ) 2 + (x x x 3 - y x y 3 - z 1 z 3 ) 2 + 
+ (x 2 x 3 - y 2 y 3 - z 2 z 3 ) 2 + {x\ + x\ + y\ + z 2 ) 2 \ 

T 4 = (yizi + y 2 z 2 + y 3 z 3 ) 2 + (x 3 y x - x x y 3 - z 2 z 4 ) 2 + (x 2 y 3 - x 3 y 2 - z x z 4 ) 2 + 
+ (xiy 2 - x 2 yi - z 3 z A ) 2 + (yf + y\ + y 2 + z|) 2 ; 



T 5 = (x 3 z 2 - x 2 z 3 - yiz 4 ) 2 + (xiz 3 - x 3 z x - y 2 z 4 ) 2 + [x 2 z\ - x x z 2 - y 3 z 4 ) 2 + 
+ (yizi + yiz 2 + y 3 z 3 ) 2 + (zf + z 2 + z 2 + zj) 2 . 

We find, in the particular case of so(5), the result from [TO], [5], that the Mishchenko integrals of even 
order (ra 2 — Ci, m 4 = K\, rue = K 2 ) are a subsystem of the Manakov integrals. 

We obtained the following list of functionally independent and Poisson commuting constants of 
motion for (|4.1[) : H, C\, C 2 ,K X , K 2 , K 3 , which proves that our system is completely integrable [2]. 

Additionally, as proved in Section 2, the system (|4.1[) admits the following "generators" integrals of 
motion: 



Fi{M) 
F 2 (M) 
F 3 (M) 
F 4 {M) 
F 5 (M) 



9 2 

4 vi 



\2 \2 1 \2 \2 ' \2 \2 1 \2 \2' 

A l ~ A 3 A l ~ A 2 A l ~ A 4 A l ~ A 5 

T2 2 ,.2 _,2 

1 1 3 1 &2 1 2 



A 2 ~ A 3 A 2 ~ A l A 2 — A 4 A 2 — A 



2 ' 



vl 



\2 \2 n \2 \2 ' \2 \2 n \2 \2' 

A 3 _ A 2 A 3 _ A l A 3 _ A 4 A 3 ~ A 5 

v\ vl vl 4 



\2 \2 1 \2 \2 1 \2 \2 1 \2 \2' 
A 4 — ^1 ^4 — A 2 A 4 — A 3 ^4 — A 5 



•/2 



\2 \2 1 \2 \2 1 \2 \2 1 \2 \2 • 

A 5 _ A l A 5 _ A 2 A 5 _ A 3 A 5 _ A 4 
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We continue by presenting the list of factorizations of the characteristic polynomials for the linearized 
equations on the tangent space to the orbit, corresponding to the equilibria from tfc n Orb Ci;C2 . 
For the equilibrium M\ b the eigenvalues are the roots of equations of the following form: 

U x t 2 + U[ = 0; Vxt 1 + V{ = 0; Wrf 4 + W[t 2 + W{' = 0; 

U x = (A 3 + A 4 ) 2 (A 3 + A 5 )(A 4 + A 5 ); U[ = 6 2 (A 3 - A 5 )(A 4 - A 5 ); 

Vi = (Ax + A 2 ) 2 (Ai + A 5 )(A 2 + A 5 ); V{ = a 2 (X 1 - A 5 )(A 2 - A 5 ); 

W! = (Ai + A 2 ) 4 (A 3 + A 4 ) 4 (A! + A 3 )(Ax + A 4 )(A 2 + A 3 )(A 2 + A 4 ); 

W{' = (Ai - AsXAj - A 4 )(A 2 - A 3 )(A 2 - A 4 )[a 2 (A 3 + A 4 ) 2 - 6 2 (A 4 + A 2 ) 2 ] 2 . 

For M 2 b the eigenvalues are the roots of equations of the following form: 

U 2 t 2 + U^ = 0; V 2 t 2 + Vj = 0; W 2 t 4 + W' 2 t 2 + W'i = 0; 

U 2 = (A 2 + A 5 ) 2 (A! + A 2 )(Aj + A 5 ); U' 2 = 6 2 (A 4 - A 2 )(Aj - A 5 ); 

V 2 = (A3 + A 4 ) 2 (A! + Aa)(Ai + A 4 ); V 2 ' = a 2 (A 4 - A 3 )(Ax - A 4 ); 

W 2 = (A 2 + A 5 ) 4 (A 3 + A 4 ) 4 (A 2 + A 3 )(A 2 + A 4 )(A 3 + A 5 )(A 4 + A 5 ); 

W!,' = (A 2 - A 3 )(A 2 - A 4 )(A 3 - A 5 )(A 4 - A 5 )[a 2 (A 2 + A 5 ) 2 - 6 2 (A 3 + A 4 ) 2 ] 2 . 

For b the eigenvalues are the roots of equations of the following form: 

f7 3 i 2 + [/^ = 0; V 3 t 2 + V^ = 0; W 3 t 4 + W' 3 t 2 + W^' = 0; 

U 3 = (Ax + A 5 ) 2 (A! + A 2 )(A 2 + A 5 ); U' 3 = -6 2 (Ax - A 2 )(A 2 - A 5 ); 

V 3 = (As + A 4 ) 2 (A 2 + A 3 )(A 2 + A 4 ); V 3 ' = a 2 (\ 2 - A 3 )(A 2 - A 4 ); 

W 3 = (Ax + A 5 ) 4 (A 3 + A 4 ) 4 (A 4 + A 3 )(Ax + A 4 )(A 3 + A 5 )(A 4 + A 5 ); 

W£ = (Ax - A 3 )(Aj - A 4 )(A 3 - A 5 )(A 4 - A 5 )[a 2 (A! + A 5 ) 2 - 6 2 (A 3 + A 4 ) 2 ] 2 . 

For M 4 b the eigenvalues are the roots of equations of the following form: 

Uit 2 + Ui = 0; Vit 2 + V( = 0: W 4 t 4 + W' A t 2 + W'l = 0; 

U 4 = (Ax + A 5 ) 2 (A 4 + A 3 )(A 3 + A 5 ); U' A = -6 2 (Ax - A 3 )(A 3 - A 5 ); 

V 4 = (A 2 + A 4 ) 2 (A 2 + A 3 )(A 3 + A 4 ); V{ = -a 2 (X 2 - A 3 )(A 3 - A 4 ); 

W 4 = (Ax + A 5 ) 4 (A 2 + A 4 ) 4 (Ax + A 2 )(Ax + A 4 )(A 2 + A 5 )(A 4 + A 5 ); 

W'l = (Ax - A 2 )(Ax - A 4 )(A 2 - A 5 )(A 4 - A 5 )[a 2 (A! + A 5 ) 2 - 6 2 (A 2 + A 4 ) 2 ] 2 . 

For b the eigenvalues are the roots of equations of the following form: 

U 5 t 2 + Ui = 0; V 5 t 2 + V^ = 0; W 5 t 4 + W£t 2 + = 0; 

U 5 = (Ax + A 5 ) 2 (A 4 + A 4 )(A 4 + A 5 ); U' b = -& 2 (A X - A 4 )(A 4 - A 5 ); 

V 5 = (A 2 + A 3 ) 2 (A 2 + A 4 )(A 3 + A 4 ); V£ = a 2 (X 2 - A 4 )(A 3 - A 4 ); 

W 5 = (Ax + A 5 ) 4 (A 2 + A 3 ) 4 (Ax + A 2 )(Ax + A 3 )(A 2 + A 5 )(A 3 + A 5 ); 

W5 = (Ax - A 2 )(Ax - A 3 )(A 2 - A 5 )(A 3 - A 5 )[ a 2 (A! + A 5 ) 2 - 6 2 (A 2 + A 3 ) 2 ] 2 . 

For M® b the eigenvalues are the roots of equations of the following form: 

U 6 t 2 + U{=0- 7 V 6 i 2 + V 6 ' = 0; W e t 4 + W^t 2 + Wg = 0: 

U 6 = (A 2 + A 4 ) 2 (A 2 + A 5 )(A 4 + A 5 ); U' t = b 2 {\ 2 - A 5 )(A 4 - A 5 ); 

V 6 = (Ai + A 3 ) 2 (Ax + A 5 )(A 3 + A 5 ); V 6 ' - a 2 ^ - A 5 )(A 3 - A 5 ); 

W 6 = (Ai + A 3 ) 4 (A 2 + A 4 ) 4 (Ai + A 2 )(Aj + A 4 )(A 2 + A 3 )(A 3 + A 4 ); 

Wq = -(A 4 - A 2 )(Aj - A 4 )(A 2 - A 3 )(A 3 - A 4 )[a 2 (A 2 + A 4 ) 2 - & 2 (A 4 + A 3 ) 2 ] 2 . 

For Mj b the eigenvalues are the roots of equations of the following form: 

U 7 t 2 + U 7 = 0: V 7 t 2 + Vj = 0: W 7 t 4 + W 7 t 2 + W 7 ' = 0; 

U 7 = (A3 + A 5 ) 2 (A! + A 3 )(Ai + A 5 ); U' 7 = b 2 (\ x - A 3 )(Aj - A 5 ); 

V 7 = (A 2 + A 4 ) 2 (A! + A 4 )(A 4 + A 2 ); Vf - a 2 (A! - A 2 )(Aj - A 4 ); 

W 7 = (A 2 + A 4 ) 4 (A 3 + A 5 ) 4 (A 2 + A 3 )(A 2 + A 5 )(A 3 + A 4 )(A 4 + A 5 ); 

W!/ = -(A 2 - A 3 )(A 2 - A 5 )(A 3 - A 4 )(A 4 - A 5 )[a 2 (A 3 + A 5 ) 2 - & 2 (A 2 + A 4 ) 2 ] 2 . 
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For b the eigenvalues are the roots of equations of the following form: 

U 8 t 2 + U[=0; V 8 t 2 + V£ = 0; W s t 4 + W^t 2 + Wg = 0; 

Us = (A 4 + A 5 ) 2 (A! + A 4 )(Ax + A 5 ); U' & = b 2 (X 1 - A 4 )(Ax - A 5 ); 

V 8 = (A 2 + A 3 ) 2 (Ai + A 2 )(Ai + A 3 ); V 8 ' = a 2 ^ - A 2 )(Aj - A3); 

W 8 = (A 2 + A 3 ) 4 (A 4 + A 5 ) 4 (A 2 + A 4 )(A 2 + A 5 )(A 3 + A 4 )(A 3 + A 5 ); 

Wg = (A 2 - A 4 )(A 2 - A 5 )(A 3 - A 4 )(A 3 - A 5 )[a 2 (A 4 + A 5 ) 2 - 6 2 (A 2 + A 3 ) 2 ] 2 . 

For b the eigenvalues are the roots of equations of the following form: 

U 9 t 2 + U{ = 0; V 9 t 2 + V^ = 0: W 9 t A + W^t 2 + = 0; 

U 9 = (Ax + A 4 ) 2 (Ai + A 5 )(A 4 + A 5 ); U' 9 = b 2 (X 1 - A 5 )(A 4 - A 5 ); 

V 9 = (A 2 + A 3 ) 2 (A 2 + A 5 )(A 3 + A 5 ); V 9 ' = a 2 (A 2 - A 5 )(A 3 - A 5 ); 

W 9 = (Aj + A 4 ) 4 (A 2 + A 3 ) 4 (A! + A 2 )(Aj + A 3 )(A 2 + A 4 )(A 3 + A 4 ); 

Wi' = (Aj - A 2 )(Ax - A 3 )(A 2 - A 4 )(A 3 - A 4 )[a 2 (Ax + A 4 ) 2 - 6 2 (A 2 + A 3 ) 2 ] 2 . 

For M*° b the eigenvalues are the roots of equations of the following form: 

U w t 2 + U[ = 0; V w t 2 + V{ = 0; W w t 4 + W{ t 2 + W{' = 0; 

U w = (A 2 + A 5 ) 2 (A 2 + A 4 )(A 4 + A 5 ); U' w = -b 2 (\ 2 - A 4 )(A 4 - A 5 ); 

V w = (Ax + A 3 ) 2 (Ax + A 4 )(A 3 + A 4 ); V{ = a 2 (Ax - A 4 )(A 3 - A 4 ); 

W w = (Ax + A 3 ) 4 (A 2 + A 5 ) 4 (A! + A 2 )(Ax + A 5 )(A 2 + A 3 )(A 3 + A 5 ); 

Wl' = -(Ax - A 2 )(Ax - A 5 )(A 2 - A 3 )(A 3 - A 5 )[a 2 (A 2 + A 5 ) 2 - b 2 (X 1 + A 3 ) 2 ] 2 . 

For the eigenvalues are the roots of equations of the following form: 

Uut 2 + U{ = 0; V n t 2 + V[ x = 0; W^t 4 + W[ x t 2 + W{[ = 0; 

U n = (A 2 + A 5 ) 2 (A 2 + A 3 )(A 3 + A 5 ); U' u = -6 2 (A 2 - A 3 )(A 3 - A 5 ); 

V u = (Ax + A 4 ) 2 (Ax + A 3 )(A 3 + A 4 ); V[ x = -a 2 (\ x - A 3 )(A 3 - A 4 ); 

W u = (Ax + A 4 ) 4 (A 2 + A 5 ) 4 (Ai + A 2 )(Ax + A 5 )(A 2 + A 4 )(A 4 + A 5 ); 

W{[ = -(Ax - A 2 )(Ax - A 5 )(A 2 - A 4 )(A 4 - A 5 )[a 2 (A 2 + A 5 ) 2 - b 2 (X 1 + A 4 ) 2 ] 2 . 

For M 12 the eigenvalues are the roots of equations of the following form: 

Ui 2 t 2 + U[ = 0; V 12 t 2 + V{ 2 = 0; W 12 t 4 + W[ 2 t 2 + W" 2 = 0: 

U 12 = (A 4 + A 5 ) 2 (A 3 + A 4 )(A 3 + A 5 ); U[ 2 = b 2 (X 3 - A 4 )(A 3 - A 5 ): 

V 12 = (Ax + A 2 ) 2 (Ax + A 3 )(A 2 + A 3 ); V{ 2 = a 2 (X 1 - A 3 )(A 2 - A 3 ); 

W 12 = (Ax + A 2 ) 4 (A 4 + A 5 ) 4 (Ax + A 4 )(Ai + A 5 )(A 2 + A 4 )(A 2 + A 5 ); 

W{' 2 = (Ax - A 4 )(Ax - A 5 )(A 2 - A 4 )(A 2 - A 5 )[a 2 (A 4 + A 5 ) 2 - 6 2 (Ax + A 2 ) 2 ] 2 . 

For the eigenvalues are the roots of equations of the following form: 

U 13 t 2 + U[ = 0; V 13 t 2 + V{ 3 = 0: W 13 t 4 + W[ 3 t 2 + W{' 3 = 0; 

U 13 = (A 3 + A 5 ) 2 (A 3 + A 4 )(A 4 + A 5 ); U[ 3 = -b 2 (X 3 - A 4 )(A 4 - A 5 ); 

V 13 = (Ax + A 2 ) 2 (Ax + A 4 )(A 2 + A 4 ); V{ 3 = a 2 (X 1 - A 4 )(A 2 - A 4 ); 

W 13 - (Ai + A 2 ) 4 (A 3 + A 5 ) 4 (Ai + A 3 )(Ai + A 5 )(A 2 + A 3 )(A 2 + A 5 ); 

W{' 3 = (Ax - A 3 )(Aj - A 5 )(A 2 - A 3 )(A 2 - A 5 )[a 2 (A 3 + A 5 ) 2 - b 2 {X l + A 2 ) 2 ] 2 . 

For M\\ the eigenvalues are the roots of equations of the following form: 

U u t 2 + U{ = 0; V 14 t 2 + V{ 4 = 0; Wut* + W{ 4 t 2 + W[' 4 = 0: 

U u = (A 3 + A 5 ) 2 (A 2 + A 3 )(A 2 + A 5 ); U{ 4 = b 2 (\ 2 - A 3 )(A 2 - A 5 ); 

V u = (Ax + A 4 ) 2 (Ai + A 2 )(A 2 + A 4 ); V{ 4 = -a 2 ^ - A 2 )(A 2 - A 4 ); 

14^14 - (Ai + A 4 ) 4 (A 3 + A 5 ) 4 (Ax + A 3 )(Ai + A 5 )(A 3 + A 4 )(A 4 + A 5 ); 

T^i" 4 - -(Ai - A 3 )(Aj - A 5 )(A 3 - A 4 )(A 4 - A 5 )[a 2 (A 3 + A 5 ) 2 - 6 2 (Ax + A 4 ) 2 ] 2 . 
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For the eigenvalues are the roots of equations of the following form: 

U 15 t 2 + U[ = 0; V 15 t 2 + V{ 5 = 0; W lb t A + W{ 5 t 2 + W{' 5 = 0; 

U 15 = (A 4 + A 5 ) 2 (A 2 + A 5 )(A 2 + A 4 ); U[ 5 = b 2 {\ 2 - A 4 )(A 2 - A 5 ); 

V 15 = (Ai + A 3 ) 2 (A! + A 2 )(A 2 + A 3 ); V{ 5 = -a 2 (X 1 - A 2 )(A 2 - A 3 ); 

W 15 = (Ai + A 3 ) 4 (A 4 + A 5 ) 4 (A! + A 4 )(Ai + A 5 )(A 3 + A 4 )(A 3 + A 5 ); 

W{> 5 = (A 4 - A 4 )(A X - A 5 )(A 3 - A 4 )(A 3 - A 5 )[a 2 (A 4 + A 5 ) 2 - b 2 {\ + A 3 ) 2 ] 2 . 

Without loss of generality, we can choose an ordering for A^'s, namely 

Ai > A 2 > A 3 > A 4 > A 5 . 



Since for the equilibria t 3 D Orb Ci;C2 , t 4 fl Orb Ci;C2 , ts n Orb Ci;C2 , tm D Orb Ci;C2 , in (~l Orb Ci;C2 , ti 3 n 
Orb Ci;C2 , ti4 fl Orb Ci;C2 and tis n Orb Ci;C2 at least one of the eigenvalues is real and strictly positive we 
have nonlinear instability for these equilibria. 

It is immediate to prove that if 

a 2 (A 2 +A 4 ) 2 ^fe 2 (A 1 +A 3 ) 2 

is equivalent from Eqs. 13.31 with 



cj [(A a + A 4 ) 4 + (Ax + A 3 ) 4 ] + c 2 [(A 2 + A 4 ) 2 + (Ai + A 3 ) 2 ] 2 , (5.1) 

and in this case we have W§ < 0. Since We > and W§ < 0, it follows that if the condition (|5.ip holds, 
then the equilibrium M® b from t 6 fl Orb Ci;C2 (together with the equilibria M^_ a b , M% _ b , M e _ a _ b ) are 
also nonlinear unstable. 

On the other hand, since 

& 2 (A 2 + A 4 ) 2 <a 2 (A 1 + A 3 ) 2 , 

(because a > b and Ai + A 3 > A 2 + A 4 ) it follows that the equilibrium M ba from t$ n Orb Ci;C2 (together 
with the equilibria M b _ a , M® b a , M^_ b _ ) are nonlinear unstable, regardless if condition (I5.ip holds or 
not. 

A similar discussion holds for the equilibria from tj Pi Orb Ci:C2 . 

We continue with the study of the remaining families of equilibria, which is tiPlOrb Ci;C2 , t 2 nOrb Ci;C2 , 
tg n Orb Ci;C2 , t 9 ("1 Orb Ci;C2 and t i2 n 0rb ci;C2 . 

For start, let us notice that energy methods using the "canonical" integrals of motion (the Hamil- 
tonian and/or the constants of motion K\, K 2 , K3) work only for a few of the above equilibria. But 
choosing particular linear combinations of some " generators" integrals of motion introduced in Section 
2 proves to be very effective in solving the stability problem with energy methods. 

To begin with, we note that the method involving linearization used above is inconclusive for the 
stability of the Mr 6 equilibrium. Therefore, we will use Arnold's method |11) . which is equivalent with 
the other energy methods [5]. 

Consider the smooth function G mn G C°°(so(5), R), where m, n are real numbers, 

G mn (M) = Fi(M) + F 5 (M) + mC x (M) + nC 2 (M). 

Choosing m, n such that dG mn (M^ b ) — 0, namely 

2b 2 2 
m - {\ 2 ~\ 2 ){a 2 -b 2 Y n -~{\ 2 -\ 2 ){a 2 -b 2 ) 

and taking into account that 

W := kerdCi(A4\ b ) n ker dC 2 (M a 1 i6 ) = span(£i, E 2 , E A , E 5 ,E 7 , E 8 , E 9 ,E W ) 
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we obtain the determinants associated to all upper-left submatrices of the Hessian d 2 G mn (M^ b )\wxW 
as follows: 



2a 2 

(A?-Al)(a*-6a) 
4a 2 [a 2 (A|-A|)+6 2 (A?-Al)] 



Dl '\2 _ \2V„2 _ J.2\ < °! 



2 (A?-Al)2(Af-Al)(a2-62) 2 

8a 2 (A 2 -A 2 )[a 2 (A 2 -A 2 ) + fr 2 (A 2 -A 2 )] 
(A?-A2)3(A?-A§)(A?-A2)(aa-&*)2 ' 
D = 16a 4 (A 2 -A 2 )(A 2 -A 2 ) 

4 ( A f-Al)4(A?-A|)(Af-Ai)(a2-6T ' 

D 32a 4 (A 2 ~A 2 )(A 2 -A 2 ) 

5 (A 2 -A 2 )5(A 2 -A 2 )(A 2 -A 2 )(a 2 -6 2 ) 2 > 
n 64a 4 (A 2 -A 2 )(A 2 -A 2 )(A 2 -A 2 ) 

6 (A? - A1)«(A? - Al)(Af - A 2 )(A 2 - %)(a* - 6 2 ) 2 > 

n = 128a 4 (A 2 -Al)(A|-Al)(Al-Al) 

7 (A 2 - A 2 )6(A 2 - A 2 )(A? - A 2 )(A 2 - A 2 )(A 2 - A 2 )(a 2 - 6 2 ) 2 < U ' 

n = 256« 4 (A 2 -A 2 )(A 2 -A 2 )(A 2 -A 2 ) 

8 (A 2 - A 2 )6(A 2 - A 2 )(A 2 - A 2 )(A 2 - A 2 )(A 2 - A 2 )(A 2 - A 2 )(a 2 - 6 2 ) 2 > U ' 

Consequently, d 2 G m „(M^ b )\wxw is negative definite, which implies nonlinear stability for the equi- 
librium M\ h . The same computations lead to nonlinear stability for the equilibria M\ _ b , M£ b and 
M!_ a _ 6 from ti n Orb Ci;C2 . 

For the remaining equilibria M ba , M_ b , M b _ a and M_ b _ a from ti PI Orb Ci;C2 a convenient energy 
function for applying Arnold's method is G mn (M) = F 4 (M) + F 5 (M) + mC\{M) + nC 2 (M). Thus, all 
eight equilibria in ti PI Orb Cl ;C2 are nonlinear stable. 

For the equilibria M 2 b , M 2 a b , M 2 _ b and M' 2 _ a _ b from t 2 P Orb Ci;C2 a convenient energy function 
for applying Arnold's method is G mn {M) = F X {M) - F A (M) +mC 1 (M) +nC 2 (M). Thus, the equilibria 
M a 2 t , M 2 a b , M 2 _ b and M 2 Q _ b are nonlinear stable. 

For the equilibria M® 6 , M^ a b , _ b and M® a _ b from t s P Orb Ci;C2 a convenient energy function 
for applying Arnold's method is G mn (M) = F 1 (M)'+ F 5 (M) + mCi(M) + nC 2 (M). For the equilibria 
M ba , M® b a , M b _ a and M® b _ a from t 8 P Orb Ci;C2 a convenient energy function for applying Arnold's 
method is G mn {M) = Fi(M) '+ F 2 (M) + mCi(M) + nC 2 (M). Thus, all eight equilibria in t 8 P Orb ci;C2 
are nonlinear stable. 

For the equilibria M b a , M® 6 a , M b _ a and M^ b _ Q from tg P Orb Ci;C2 a convenient energy function 
for applying Arnold's method is G mn {M) = F A (M) - F 5 (M) + mC 1 {M) +nC 2 (M). Thus, the equilibria 
M b a , M® b a , M b _ a and M® 6 _ a from t 9 P Orb ci:C2 are nonlinear stable. 

For the equilibria M* 2 , M^ 2 b , M^ 2 _ b and A/1 2 _ b from ti 2 P Orb ci;C2 a convenient energy function 
for applying Arnold's method is G mn (M) = F 1 (M) - F 3 (M) + mC x (M) + nC 2 {M). For the equilibria 
M 12 , M^ 2 , M b 12 _ a and M^ 2 _ from t 12 P Orb ci;C2 a convenient energy function for applying Arnold's 
method is G mn (M) = F 3 (M)'- F 5 (M) + mC 1 {M) + nC 2 (M). Thus, all eight equilibria in t 12 P Orb ci;C2 
are nonlinear stable. 

For the equilibria M b a , M 2 b a , M b _ a and M 2 b _ a from t 2 P Orb ci;C2 and for the equilibria M® b , 
Mt a b , _ b and M^_ a _ b from tg P Orb Ci;C2 the stability problem remains open and it is likely that 
a bifurcation phenomenon occurs. Such a phenomenon appears for some equilibria in the case of the 
so(4) free rigid body and is extensively studied in [6]. 

We proved the following result: 
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Theorem 5.1. (i) The equilibria from t 3 n Orb Ci;C2 , t 4 fl Orb Ci;C2; i$ n Orb Ci;C2 , tio H Orb Ci;C2 , tn n 
0rb Ci;C2 , ti3 n Orb Ci:C2 , in n Orb Ci;C2 and ti5 n Orb Ci;C2 are unstable. 

(ii) T/ie equilibria M® a , M fc 6 _ a; M^ 6 a and M® b _ a from t 6 n Orb 

ci;c 2 anc ^ respectively the equilibria 
M b,a> M b-a> M -b.a and M -b,-a from *7 n 0rb c i;C2 are unstable. 

(iii) a) If the condition 

c\ [(Ax + A 3 ) 4 + (A 2 + A 4 ) 4 ] ^ c 2 [(Ai + A 3 ) 2 + (A 2 + A 4 ) 2 ] 2 

holds, then the equilibria M® b , M 6 _ a b , M® _ b and M^_ a _ b from i§ fl Orb Ci:C2 are unstable, 
b) If the condition 

c\ [(A 2 + A 4 ) 4 + (A 3 + A 5 ) 4 ] ^ c 2 [(A 2 + A 4 ) 2 + (A 3 + A 5 ) 2 ] 2 
holds, then the equilibria M' J ab , M 7 _ a b , M 7 a _ b and M 7 _ a _ b from n Orb Ci;C2 are unstable. 

(iv) The equilibria from ti fl Orb Ci;C2 , tg fl Orb Ci;C2 and ti 2 (~l Orb Ci;C2 are nonlinear stable. 

(v) TTie equilibria M 2 b; M 2 a 6 , M 2 _ fc and M 2 a _ b /rom t 2 n Orb 

ci;c 2 anc ^ respectively the equilibria 
Mj a , Q , M b _ a and M 9 _ b _ a from tg n Orb Ci;C2 are nonlinear stable. 

Acknowledgements. The author acknowledges the very helpful discussions with Professors Tudor 
Ratjiu and Petre Birtea. 

References 

[1] A. S. Mishchenko, Integral geodesies of a flow on Lie groups, Fund. Anal. Appl. 4 (1970) 232-235. 

[2] S. V. Manakov, Note on the integration of Euler's equations of the dynamics of an n-dimensional 
rigid body, Fund. Anal. Appl. 10 (1976) 328-329. 

[3] T. S. Ratiu, The motion of the free n-dimensional rigid body, Indiana Univ. Math. J. 29 (1980) 
609-629. 

[4] L. Feher and I. Marshall, Stability analysis of some integrable equations for SO(n), J. Nonlin. 
Math. Phys. 10 (2003) 304-317. 

[5] A. Spiegler, Stability of Generic Equilibria of the 2n-Dimensional Free Rigid Body Using the 
Energy-Casimir Method, thesis, University of Arizona, 2004. 

[6] P. Birtea, I. Casu, T. Ra^iu, M. Turhan, Stability of equilibria for the so (4) free rigid body, 
arxiv.org/abs/0812.3415. 

[7] A. V. Bolsinov and A. T. Fomenko, Integrable Hamiltonian Systems (Chapmann & Hall/CRC Press, 
2004). 

[8] P. Birtea and M. Puta, Equivalence of energy methods in stability theory, J. of Math. Phys. 48 
(2007) 042704. 

[9] C. Morosi and L. Pizzocchero, On the Euler equation: Bi-Hamiltonian structure and integrals in 
involution, Lett, in Math. Phys. 37 (1996) 117-135. 

[10] M. Ugaglia, Sistemi dinamici integrabili su algebre di Lie: funzioni di Casimir e significato del 
paramctro spettralc, thesis, University of Torino, 1994. 

[11] V. I. Arnold, Conditions for nonlinear stability of stationary plane curvilinear flows of an ideal 
fluid, Doklady 5 (1965) 773-777. 



14 



